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ABSTRACT 


The diffusion of protons in the outer radiation belt due to viola- 
tion of the third adiabatic invariant has been examined. The particular 
mechanism studied is one where variations in the solar-wind intensity 
produce magnetic disturbances causing motion of particles between 
L -shells. A Fokker-Planck diffusion equation is used with terms 
describing Coulomb energy degradation and charge -exchange losses 
with the ambient atmosphere. A source near the magnetopause is as- 
sumed. The diffusion constant is numerically evaluated by an analysis 
of observed frequencies of sudden impulses and commencements. The 
equation is solved to obtain radial and energy distributions. A com- 
parison of the results with the measured proton fluxes of Davis, Hoff- 
man, and Williamson indicates that a diffusion process of this general 
nature seems to be responsible for the observed outer-belt protons. 
For the parameters selected, the magnetic disturbances studied may 
not be sufficient to produce the required diffusion rate. The effects of 
possible changes in parameters are discussed. 



Introduction 


Recent observations by Davis and Williamson [l963] and Davis 
et al. [l9 64] in the L = 2 to 8 earth radii region of the outer radiation 
belt confirm the existence of large fluxes of 0.1 to 10 Mev protons. 
Except for the more energetic protons beyond L = 5, the intensities 
appear to be very stable over times of the order of years. An out- 
standing feature of these protons is the large but smooth variation in 
their spectra with L and equatorial pitch angle. 

The recent theoretical study of the spectra of these protons by 
Dungey et al. [l9&5] strongly suggests that the source is near the edge 
of the magnetosphere and that some process that violates the third 
adiabatic invariant for trapped particles (i.e., permits motion between 
L-shells) plays an important role in populating the radiation belt. 

Parker [i 960 ], Davis and Chang [ 1962 ], and Tverskoy [1964] have 
studied one such process which must surely operate according to present 
interpretations of geomagnetic fluctuations due to solar plasma. 

This mechanism operates in the following way. Consider particles 
trapped near the equator, which drift around the earth on a constant 
magnetic field path. When the solar wind intensity is low and the geo- 
magnetic field is relatively undisturbed, their drift paths may be 
represented by the solid line of Figure 1(a). If the solar wind intensity 
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increases in times short compared to the longitudinal drift period of 
the particles, the particles will follow the lines of force as they are 
compressed during the magnetic disturbance. They will move into a 
smaller ring which is displaced away from the sun, as indicated by the 
dashed line in Figure 1(a) and the solid line in Figure 1(b). Different 
portions of the displaced ring will follow different constant magnetic 
field paths (dashed lines in Figure 1(b)) during subsequent longitudinal 
drift. In general, these paths will be closer to the earth on the dark 
side. Now if the solar wind decays to its original low value in times 
long compared to the longitudinal drift period, the third adiabatic 
invariant will remain valid for the particles during this slow expansion 
phase. The net effect is that any initially narrow ring of particles 
(dashed line of Figure 1(d)) will become broadened (between solid lines 
in Figure 1(d)) and diffusion in L- space will occur. . 

Sudden impulses and sudden commencements have rise times of 
a few minutes and decay times of a few hours. Typical protons seen by 
Davis et al. [l964] have drift times of the order of a half hour or so. 

(The drift time of equatorially -trapped, non- relativistic particles is 
44/LE minutes [ Lew , 1 9 6 1] , with L in earth radii and E in Mev.) Thus 
the protons considered here should be diffused by this mechanism. 

Parker [i 960 ] evaluated the mean square broadening per magnetic 
disturbance and obtained the large -L portion of the steady state distribution 
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Figure 1 -Schematic drawing of a sudden c< 
during corresponding phases of the sud< 



of electrons by using a diffusion equation and considering many dis- 
turbances. Davis and Chang [l 9 6 2 ] found that for this mechanism the 
Fokker-Plank equation was appropriate |see also Parker [ 1963 ], Tverskoy 
[l964], and Dungey [l 9 5 ] ^ . They concluded that the mean displacement 
of a ring of particles, <(Ar> , had a large effect on the equatorial distri- 
bution. Davis and Chang [l 9 6 2 ] also obtained steady -state radial dis- 
tributions by assuming a source near the magnetopause and assuming 
that the particle density vanished at a particular value of L near the 
earth. This latter assumption is quite reasonable, since the diffusion 
rate becomes small very rapidly as the earth is approached. However, 
unless specific loss rates are considered, it is not possible to de- 
termine where the density vanishes and, thus, the effectiveness of the 
mechanism. 

Tverskoy [1964] added Coulomb energy loss to the Fokker-Plank 
diffusion equation, obtained solutions for relativistic electrons, and 
concluded that this mechanism was important for these electrons. He 
also compared diffusion times with Coulomb energy-loss lifetimes of 
protons and concluded that if sources existed, Mev protons and a par- 
ticles should exist in the outer belt. We also made a similar compari- 
son [ Mead and Nakada , 1964]. Although this procedure does give an 
indication of the effectiveness of the mechanism, it is rather unsatis- 
factory, since the interpretation of "diffusions times" and Coulomb 


4 



energy-loss lifetimes is vague and dependent on the proton 
spectrum. 

In the present paper, therefore, we have obtained solutions to 
the Fokker-Planck equation containing terms describing energy losses 
due to Coulomb interaction with the ambient electrons as well as charge 
exchanges losses with the neutral atmosphere. Both radial and energy 
distributions are obtained. 

Instead of an image dipole model of the disturbed field, we used 
the more accurate magnetosphere model of Mead [l964], based on the 
solution to the Chapman-Ferraro problem of a perfectly -conducting 
plasma incident upon a three-dimensional dipole field. Also, sudden 
commencements and sudden impulses were interpreted in terms of the 
compression of the magnetosphere from an initial quiet-time con- 
figuration, caused by an increase in the solar wind intensity, rather 
than assuming that the solar wind was absent between magnetic dis- 
turbances, as did Parker and Davis and Chang. To obtain a numerical 
value for the diffusion coefficient in the Fokker-Planck equation, we 
used an estimate of the actual number of observed sudden impulses and 
commencements of various amplitudes during the period 1958 to 1961. 

The solutions to the equation were then converted to integral 
fluxes j ( >E ) corresponding to the seven energy thresholds of Davis 
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and Williamson's detector. The resulting theoretical curves were 
compared with their experimental fluxes. The curves are remarkably 
similar, indicating that a diffusion process of this general nature seems 
to be responsible for the observed outer belt proton fluxes. Finally, 
the effects of modifying the various coefficients in the Fokker-Planck 
equation are discussed. 
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Steady -State Solutions 


Many of the simplifying assumptions that have characterized 
previous calculations have been retained. The direction of the solar 
wind is assumed to be perpendicular to a dipole geomagnetic field. 

Only near -equatorial protons are considered. The magnetic moment 
and longitudinal invariants for the protons are assumed to remain 
constant during L-motion. Pitch angle changes due to the ambient at- 
mosphere, which should be small, have been neglected. A steady 
source has been assumed near the quiet time magnetopause, which we 
assume is at 10 earth radii in the solar direction. We have also as- 
sumed that the mechanism applies out to distances quite near the 
magnetopause. 

We consider a distribution function, y(/x, r, t) d/xdr, equal to the 
number of protons in dr at geocentric radius r, with magnetic moments 
in the interval d/x at fi, and with equatorial pitch angles (EPA) between 
77/2 - § 0 (r) and 77/2 + § 0 (r) , where S Q (r) is a small angle whose changes 
with r are determined by the preservation of the magnetic moment and 
longitudinal invariants. This distribution function is essentially equiva- 
lent to Davis and Chang's <fi* . The magnetic moment /x - E/B = r 3 E/M, 
with energy E in Mev has been used. With r in earth radii, M = 0.312 
gauss. In Figure 2, p. is plotted versus E for a number of r values. 
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The magnetic moment, /x, is used in the distribution function, rather 
than energy, because as particles move in L-space, /x is constant except 
for Coulomb losses, whereas the energy varies with r. 

The Fokker-Plank equation with loss terms is then [Davis and 
Chang , 1962; Parker , 1963]: 



The angular brackets denote averages over many storms and r cx is the 
1/e charge -exchange lifetime. For steady-state solutions, dy/dt = 0. 

In the evaluation of the averages in the first two terms on the 
righthand side of equation (1), ^Ar^ and <((Ar) 2 ^> have first been eval- 
uated for a single sudden commencement; then, time averages have 
been obtained by considering disturbances of all sizes and the rate of 
occurrence of these disturbances. For any single event, the deforma- 
tion of the magnetic field by the solar wind is used to evaluate <(Ar^> and 
<^(Ar) 2 ^. Mead, [l 964] and Midgley [l964] have described this deforma- 
tion with spherical harmonic expansions. They find that the distorted 
field is adequately described out to distances quite near the magnetopause 
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by the use of two coefficients in the expansion. Higher-order terms in 
the more precise description of the distorted field have been found to 
have little effect on this mechanism. 

Since the image dipole description of the distorted field as used 
by Parker [i 960 ] and Davis and Chang [l 962] is a special case of the 
two-coefficient spherical harmonic expansion, their results for <^Ar> and 
^(Ar) 2 ^> have been used, with changes in only the values of the harmonic 
coefficients. 

In prior calculations of <^(Ar) 2 ^>, the solar wind was assumed to 
be absent between world-wide magnetic disturbances. Recent measure- 
ments of both solar wind parameters [Snyder and Neugebauer , 1964] and 
the magnetospheric boundary [ Ness et al., 1964] indicate that even 
during undisturbed times, a boundary is present at about 10 earth radii 
in the solar direction. In accordance with these findings, the quiet -time 
boundary, r bo , is assumed to be at 10 earth radii. 

In terms of the first two constants in the spherical harmonic 
description of the distorted field, the equatorial field is given by 


B 



r cos 4 > 


( 2 ) 
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with similar, but more complicated expressions for the distorted field 
off the equator. Here r and r b (the distance to the boundary in the 


solar direction) are in earth radii, M = 0.312 gauss, and <p is the longi- 
tude measured from the local noon meridian. For the image dipole 
model used by Parker and by Davis and Chang, a l - M/8 and a 2 = 3M/16 
In Mead's more accurate magnetosphere model, a x = 0.81 6M - 0.2515 
gauss and a 2 = 0.673M = 0.210 gauss ^compare equations (10-12) of 
Mead [l 9 64] with equation (6) of Davis and Chang [ 1962 ]^. 


An analysis similar to Parker's and Davis and Chang's gives 


<W> = 1 (t) 

which reduces to previous results if r bo is very large. Since for a 
worldwide disturbance of average equatorial size AB, 



AB = 




) 
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the mean square displacement may be written 


<(Ar ) 2 > = 



( 4 ) 


Because of the effect of induced currents in the earth and ionos- 
phere, the disturbance in the horizontal component which is actually 
observed, AB' , is about 50% larger than that produced by the compression 
of the magnetosphere, AB. Listed in Table 1 are a compilation through 
the courtesy of Dr. M. Sugiura of the average number of sudden impulses 
and sudden commencements per year of observed size AB' for the years 

Table I 

Observed frequencies of sudden commencements and sudden 


impulses of various sizes during the period 1958-1961. 


AB' (gammas) 

Number/ Year 

Relative Effect 

> 100 

0.5 

.20 

60 - 100 

1.8 

.25 

40 - 60 

2.3 

.10 

20 - 40 

21 

.29 

5 - 20 

61 

.12 

2 

720 

.04 
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1958 to 1961. (A sudden commencement is followed by the main phase 


of a magnetic storm; a sudden impulse is not.) 


Using the values of Table I and setting AB = 2/3 AB' , and using 
Mead's coefficients for and a 2 , the diffusion coefficient <^Ar)^At)> 
has been evaluated with Equation (4) to give 



031r -fe) 


10 


(earth radii ) 2 / /day 


(5) 


Davis and Chang [l962j have derived a relationship between (Ar)> and 
<(Ar) 2 )> : 


<Ar> _ g_ 
<(Ar) 2 > " 2r 


( 6 ) 


where g is a constant that appears to be independent of the harmonic 
coefficients. They find g = 8; when his equations are recast in similar 
form, Tverskoy [l 9 64] finds g = 6.5. In the present calculations, a 


13 


range of g has been used. Equation (6) has been used to obtain: 



.031 

2 



earth radii/day. 


( 7 ) 


For the third term on the right in Equation (1), the Coulomb 
energy loss formula has been used: 


_ 4tt e 4 p 'CnA 

<9t m v 


where p is the electron density, / CnA = 20, m is the electron mass, and 
v is the proton velocity. Thus 



dp <?E 
<9E dt 


- 477 -f.nA 


e 4 p 

m 




= - 3. 55 x 10" 6 p 


l/? Mey 

y p gauss -day 


( 8 ) 
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with r in earth radii and p in electrons / cm 3 . For the electron density 


we have used 


P 


(8 x 10 3 \ _ 3 

+ 50 ) cm 


For the charge -exchange term, 


1 

r 

cx 


VP 0 V 


where p Q is the neutral hydrogen density, which we have taken to be 


7. 35 x 10 3 


cm 


with r in earth radii. For the charge -exchange cross-section, a, 
values given by Bates and McCarroll [1962] were used. The choice 
of electron and neutral hydrogen densities is treated below in the 
Discussion section. 
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When these terms are put into equation (1) and x = r / r bo> Equation 


(1) becomes: 


+ ■ ^ - 0 (9) 

which is the same as Equation (30) of Davis and Chang [l 9 62 ] , except 
for the loss terms. 

This equation has been solved numerically on a computer with 
boundary values at small x , larger, and at x = 1. At small x , which 
corresponded to locations within the dense atmosphere of the earth, 
dy/dx = 0 has been used for all fi . This condition is equivalent to 
y = 0 close to the earth and gave results that were independent of the 
location of the inner boundary as long as it was well within the dense 
atmosphere. 

At large /x, which corresponds to the higher energies at which 
protons have negligible charge exchange loss, the boundary conditions 
were obtained by solving Equation (9) without the charge exchange 
term. Except near the /x at which this boundary condition was applied, 
the solutions obtained were found to be insensitive to the boundary 
values that were used. 


031 


dx 1 


(* 10 y) _ B U 9 y) 
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At x = 1, which corresponds to the magnetopause and the as- 
sumed source location, a source energy spectrum was used. Since the 
/x spectrum was found to be essentially independent of r except near 
the earth, a spectrum similar to those measured by Davis et al. [l 9 64] 
near 4 earth radii was used. This spectrum is of the form e“ M/ ^o. The 
qualitative results for different choices of /x 0 were similar; however, 
the quantitative results were strongly dependent on the source spectrum. 

Figure 3 shows a distribution function plotted versus r for /x = 200 
Mev/gauss, for which the proton energy is 1 Mev at 4 earth radii. 
Coulomb loss is the only important loss process at this /x value. The 
dashed curve in Figure 3 was obtained by Davis and Chang [1962] with 
the same g = 8 but by arbitrarily setting y = 0 at r = 3. Since the 
diffusion parts of the differential equations are the same, the two 
curves are similar at large r, where diffusion should be dominant, and 
differ at small r where losses become important. Their loss rate is 
zero for r > 3 and is infinite at r = 3, so their curve drops very rapidly 
near r = 3. The loss rate used in these calculations for \x = 200 is 
almost constant for all r values; the effect of this distributed source 
and sink is apparent in the flattened curve with slower decrease at 
small r . 
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DAVIS AND 
)T CHANG 


r (earth radii) 


Figure 3-Comparison of the distribution functions 
of Davis and Chang ( 1962 ) and corresponding cal- 
culated results of this paper at ji~ 200Mev/gauss, 
with g = 8- The curves nave been normalized to 
each other at r =10 earth radii. At large r, where 
diffusion is dominant, the curves are very similar; 
at smaller r, differences are due to the different 
treatment of loss processes. 




In Figure 4, y is plotted versus r for a number of /x values. For 
fx greater than about 100, the curves are very similar in shape. The 
smaller /x curves peak at slightly larger r and drop more rapidly at 
small r. For p. = 10, the dip at r = 3.5 is due to charge -exchange loss; 
for r < 3.5, there is an additional, noticeable supply of protons from 
larger fj. values due to energy degradation. 

Figure 5 shows the gradual change in the /j. spectrum as protons 
diffuse inward. Except at small /x, the r = 6 curve has the same ex- 
ponential form as the assumed source spectrum at r = 10, shown by the 
dashed line. At r values less than 4, the /x spectra show a gradual 
hardening with smaller r. 

In Figures 4 and 5, y has been plotted with fj. as a variable because 
this choice best illustrates the trends in the solutions. In the Appendix 
y/r is shown to be proportional to the differential flux, j (protons /cm 2 - 
sec -ster-Mev). In Figure 6, relative differential fluxes are plotted 
versus r for various proton energies. In Figure 7, energy spectra are 
given at different r . Figures 6 and 7 do not have the same overall 
normalization as Figures 4 and 5. 

The results shown in Figures 4-7 are typical in that variations 
in parameters such as g, the source energy spectrum, diffusion rate, 
and the atmosphere gave qualitatively similar results. The effects of 
varying the parameters are discussed in the next two sections. 
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r (earth radii) 


Figure 6— Relative fluxes versus geocentric distances for various proton energ 
curves are transformed from those in Figure 4. 
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Comparison with Measurements 


Experimental results on proton fluxes in the outer belt have been 
obtained with threshold detectors [Davis and Williamson . 1963] which 
gave fluxes, j(>E), in number of protons per cm 2 -sec-ster with energy 
greater than E. The corresponding fluxes have been calculated from 
our theoretical curves and are shown to the right of Figure 8. These 
results have been normalized to the peak flux from measurements. To 
the left in Figure 8 are the fluxes for near -equatorial particles that 
have been measured by Davis et al. [1964]. 

This comparison shows good qualitative similarity between the 
satellite measurements and the results of this model. One important 
quantitative difference is that the calculated curves are displaced out- 
ward by about 0.7 earth radii. This means that for the parameters 
that have been used, this mechanism is not effective enough to explain 
the measured fluxes. 

Because parameters may be re-evaluated, a study has been 
made to see what changes in them would be required to shift the cal- 
culated curves toward lower r . In this study, the atmosphere has been 
kept constant, and the source spectrum has been adjusted to keep the 
ratio of peak fluxes for various threshold energies in agreement with 
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measurements. The diffusion rate, which is proportional to the constant 
in Equation (5) and the same constant in the first term of Equation (9), 
was increased for various values of g. Table II gives the factor of 
increase in the diffusion rate that is necessary to shift the calculated 
results enough to bring about agreement in r with the measured curves. 

Table II 


Increases in the diffusion constant in the first term of Equation (9) 
which are needed to give good agreement with the measured 

fluxes of Figure 8 


g 

Increase in 
Diffusion rate 

4 

x 2 

5 

x 3.2 

6 

x 5. 5 

8 

x 8 
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Discussion 


The results shown in Table II indicate the necessity for an in- 
crease in the diffusion rate over the value that was calculated in Equa- 
tion (5) unless g is small. In this section, the various factors that 
entered into the evaluation of the diffusion rate are examined, and other 
parameters such as g and the atmosphere are discussed. 

Let us assume that the numbers and magnitudes of clearly - 
identifiable sudden commencements and sudden impulses as given in 
Table I are correct. The observed magnitudes at the earth's surface 
are likely to be lower-limit values, since there may be attenuation of 
magnetic disturbances on their passage through the ionosphere. A 
comparison between satellite and ground measurements, such as has 
been made by Nishida and Cahill [l 9 64] , gives some indication of this, 
but more measurements are necessary. Measurements within the 
radiation belt of these sudden changes would be ideal for evaluating 
both numbers and magnitudes. Since the diffusion rate varies ap- 
proximately as (AB ) 873 , an attenuation of 30% would give an increase 
in the diffusion rate by about a factor of 2. 

Table I includes only positive sudden commencements and impulses. 
Nishida and Cahill [l964] state that negative sudden impulses occur about 
half as often as positive sudden impulses. Since Parker's [i960] 
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analysis indicates that negative impulses are as effective as positive 
ones in producing diffusion, the inclusion of the negative impulses would 
increase the diffusion rate by about 1.5 if the size distribution of nega- 
tive impulses is the same as positive ones. 

Table I includes only clearly -identifiable sudden impulses and 
commencements. But any world-wide sudden change in times of less 
than about 5 minutes that is followed by an unchanged field or slow 
decay times of about 20 minutes contributes to the diffusion rate. An 
examination of satellite records may be necessary to evaluate the 
contribution of such sudden changes that are not identifiable at the 
ground as sudden impulses or commencements. Some of these events 
are likely to produce energy and r dependent changes in the diffusion 
rate. The contribution of these events to the enhancement of the dif- 
fusion rate is at present unknown. 

Falthammar [1965] has developed a method for evaluating the 
effects of very general magnetic disturbances on the diffusion rate, 
whereby one makes a power spectrum analysis of the disturbance. 

Using this method, it would be possible to include contributions due to 
magnetic disturbances other than sudden commencements and impulses. 
Such an analysis could conceivably increase the diffusion rate sub- 
stantially. Falthammar mainly considered the possibility that time 
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variations in electric fields would contribute to diffusion in the outer 


radiation belt. Unfortunately, measurements of such electric fields and 
their variations do not at present exist. 

In the calculation of the diffusion rate, the quiet-time boundary, 
r bo , was assumed to be at 10 earth radii. Changes in the assumed 
boundary have a relatively small effect on the diffusion rate. An in- 
crease of 2 earth radii in the quiet-time boundary location reduces the 
diffusion rate by about 20 percent, while a decrease of 2 earth radii 
increases the diffusion rate by about the same amount. 

The values of the two coefficients in the spherical harmonic 
analysis of the distorted magnetic field have a relatively large effect 
on the diffusion rate. As an example, if the image dipole coefficients 
had been used, the diffusion rate would be about eleven times more 
rapid than the rate for the coefficients that were used. The coefficients 
used give good agreement with the measurements of the shape of the 
magnetopause [ Ness et al., 1964; Mead and Beard , 1964]. However, 
they are deficient in explaining the noon-midnight shift of trapped 
energetic electrons that are observed by polar -orbiting satellites at 
low altitudes and high latitudes. Wiliams and Mead [l965] find that 
good agreement with measurements can be obtained if, in addition to 
the magnetopause current system, a current sheet which produces a 
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magnetospheric tail field is added. In the closed trapping region of the 
magnetopause, the effect of this current sheet is roughly equivalent to 
increasing the second harmonic coefficient. About a factor of two in- 
crease over values used here are necessary to give field distortions 
that are in agreement with the trapped electron asymmetry. If the 
current system that produces the magnetospheric tail also responds 
to sudden changes in the solar wind, then a factor of 2 increase in a 2 
would give a factor of 4 increase in the diffusion rate, since a 2 appears 
squared in Equation (4). There is evidence for increases in the tail 
field during periods of magnetic activity [Ness , private communication]. 
This response of the tail field to the solar wind appears plausible, 
since the tail field is confined by the solar wind, and pressures within 
the tail depend upon solar wind parameters. 

The parameter g is likely to be near the value of 6.5 or 8 that 
Tverskoy [l964] and Davis and Chang [1962] obtained. However, a 
range of g values have been used in these calculations since: (1) in 
determining g, they only considered a transition from no solar wind 
to a strong solar wind instead of from a finite solar wind to a stronger 
solar wind; (2) they did not consider negative sudden events or the 
possible coupling of these negative events with positive events; and (3) 
the results of the two studies are different. The g parameter was 
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found to have a large effect on the /i. distributions as in Figure 3 and 4, 
but a much smaller effect on the shape of integral flux curves like those 
in Figure 8. Thus, it does not appear possible to determine g from a 
comparison of calculations with experiment. 

As is to be expected, both the electron and neutral hydrogen 
densities have a large effect on calculated results. If the densities 
are appreciably larger then those assumed here, then even larger in- 
creases in the diffusion rate will be necessary if the mechanism under 
consideration is to contribute appreciably to the explanation of the outer 
belt protons. 

For Coulomb energy loss, only the electron density above 1.5 
earth radii was found to be important. The profile that was used was 
obtained by using an average density at 1000 km altitude as obtained by 
Alouette [Thomas and Sader , 1964; Bauer and Blumle , 1964] of 6000/cm 3 , 
and an average density of 500/cm 3 at 2 earth radii that is due to Bowles 
[1963]. An average ion mass of 6 at 1000 km was assumed with a de- 
crease to 1 at 2 earth radii; a temperature profile due to Serbu [l 9 65] was 
also assumed. A shift in the density profile to larger or smaller altitudes 
gave a corresponding shift in the calculated curves of Figure 8. 

The neutral hydrogen density has its only appreciable effect on the 
lower energy protons, for which the charge -exchange cross-section is 
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large. The four lowest-energy-threshold curves of Figure 8 were ap- 
preciably altered by changes in the high-altitude density profile. Changes 
in the density profile below about 4 earth radii, however, had little ef- 
fect on the results. The profile that was used corresponds to the 1000°K 
one due to Thomas [l 9 6 3] with an exospheric base density of 10 4 /cm 3 
at 520 km altitude. This density was reduced by a factor of 2 beyond 
r = 2, in accordance with recent calculations by Liwshitz and Singer 
[1965]. 
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Source Requirements 


Because source requirements may help determine whether this 
mechanism plays an important role in populating the radiation belts, 
they have been estimated. These estimates have been made with the 
assumption that the diffusion rate is rapid enough to give agreement in 
the radial position of the peak fluxes between calculated and experi- 
mental results. The variation of estimated source requirements with 
changes in the diffusion rate and with g, however, were relatively small. 

Near the subsolar magnetopause, an energy spectrum for the 
-e/e* 

source of the form e , with E Q in the 7 to 10 kev range, gives fairly 
good agreement with experiment. A 2 to 3 percent addition of an ex- 
ponential spectrum with E 0 near 90 kev improves the agreement of the 
calculated 1.688 Mev threshold fluxes with experiment. A power law 
spectrum of the form E" a , with a near 5, also gives relative fluxes 
that are in good agreement with experiment. 

The source strength requirement was estimated by calculating 
proton loss rates due to charge exchange and energy loss. First the 
relative total number of protons in dr at each r was calculated by 
evaluating Jy d^ over ail fi. The results are shown in the top curve of 
Figure 9. Then the relative charge -exchange loss rate per day at each 
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r (earth radii) 

in the outer belt. The top curve Jy (fi, r, t) d/x # 
f protons in dr as a function of r. The lower solid 
ge-exchange loss rate per day. The lower dashed 
ill fi, gives the Coulomb loss rate per day. 




r was evaluated by calculating J ( y/ T cx ) d/x. The lower solid line of 
Figure 9 shows these results. 

In general. Coulomb energy loss acts as a source for lower (jl 
values. But for very small /x the degradation of energy is so rapid that 
the distribution function goes to zero rapidly. The product - <^x/At)>y, 
however, remains essentially constant at small /x, and is equal to the 
Coulomb loss rate per day at each r. This product was evaluated at 
small fjL 9 and was found to be relatively insensitive to the /x value chosen. 
The results are shown as a dashed line in Figure 9. 

From these curves, it is evident that most of the loss occurs 
due to Coulomb energy degradation at large r . By comparing the 
integral over r of the relative loss rates with the integral over r of the 
total relative number of protons, we get the result that 1.2 percent of 
the trapped protons are lost per day. The corresponding loss rate for 
r <6 only is about an order of magnitude less. 

Assuming that the solar wind is the source of these protons, the 
fraction of solar wind protons which must be supplied to equal this loss 
rate has been estimated. The protons were assumed to fill a torus with 
an internal radius of 4 earth radii and a radius of 6 earth radii between 
the center of the earth and the locus of the internal radius. Proton 
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fluxes were estimated from experimental results for small pitch angles. 
The solar wind injection area was assumed to be a circle with a radius 
of 4 earth radii, and a solar wind flux of 10 8 protons/cm 2 -sec was 
also assumed. With these numbers, it was found that on the average 
approximately 10" 6 of the solar wind protons must be supplied to the 
magnetosphere in order to replace those lost due to charge exchange and 
energy loss. 

Whether or not sufficient protons of the required energies exist 
in the solar wind has not been determined. If the measured tempera- 
tures of about 10 ev and directed energies of about 1 kev in the inter- 
planetary region are assumed, then unaltered solar wind protons do not 
have the necessary energies and spectrum to supply source require- 
ments. It may be that the necessary changes occur in the transition 
region between the shock front and the magnetopause. 

If more energetic protons from the sun or interplanetary space 
are the source particles, time average fluxes exceeding 10 2 /cm 2 -sec 
are required. 
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Conclusions 


These studies strengthen the suggestion by Dungey et al. [l 965 ] 
and Tverskoy [l 964] that the source of the outer belt protons is near 
the magnetopause and that some process that violates the third adiabatic 
invariant is important in moving protons into the region of the outer 
radiation belt. Although uniqueness cannot be claimed, the similarity 
of the calculated fluxes to measurements gives support to the idea that 
some diffusion process is important. 

Conclusions on whether or not this diffusion is caused exclusively 
or in major part by magnetic disturbances due to solar wind fluctuations 
must await further experimental and theoretical results. Similar con- 
clusions may be drawn about the source and injection mechanism. 

Frank et al. [l 9 64] , Frank and Van Allen [1964] and F rank [l 9 6 5] 
have reached similar conclusions on the importance of some diffusion 
process in their studies of energetic electrons (E > 1.6 Mev) in the 
outer radiation belt. They further find that large variations in the 
source of electrons are likely. In contrast, the source for protons is 
likely to exhibit considerably smaller time variations, since proton 
fluxes at r < 5 earth radii are quite stable. 
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Although proton fluxes have been quite stable over the last few 
years, changes are likely to occur during other phases of the solar 
cycle if the processes considered here are important. Changes in 
parameters such as source strength, source spectrum, magnetic 
activity, neutral hydrogen density, and the electron density can cause 
predictable changes in proton fluxes and spectra. It should be inter- 
esting to follow changes in measurable parameters and proton fluxes 
as solar maximum is approached. 
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Appendix 


In this appendix, the relationship between the distribution function 
y( fjL, r) d/xdr and the directional flux, j(E, r, a) dEdAdftdt is derived. 
Davis and Chang [ 1962 ] have obtained a similar expression but with 
slightly different assumptions about the nature of the fluxes. Here j at 
the equator is assumed constant at any r over the pitch angle range 
from 77/2 ~ § 0 to n /2 + S Q where S Q is a small angle whose change with 
r is determined by the preservation of the magnetic moment and the 
longitudinal invariant. Only those particles with equatorial pitch angles 
between 77/2 - S Q and 77/2 + § 0 contribute to y. 

Thus y is given by: 


r +x m r 1T '' 2+s j dE 

y - 277r - j- 277 sin 6 d<9 r dk (1A) 

J -\ a Jv/2-S 


where 0 is the pitch angle, S(X., r) varies with K the latitude, since 
pitch angles change with k, and ±^ m corresponds to the mirror latitude 

i 

for a proton with equatorial pitch angle, 77/2 ± § 0 . 
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The integration over 8 gives the result: 


y 



2 sin S dk 


From the magnetic moment invariant, the following are obtained: 


sin2 (i ~ s ) 

sir,2 (i - § o) 




Vl + 3 sin 2 k 

f 

cos 6 k 

m 



for X small. Also 

m 


sinS = {K 2 ~ k2 ) l/2 


for X and X small. 

m 


(2A) 
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Using these results, Equation (2A) becomes: 


* 


y 


C r 2 


i <kL 

v dE 


X 2 

m 


(3A) 


where C is a constant. 


When X m is small, the longitudinal invariant J is: 


J 


4 



r dA. 


= K v r X 


2 

m 


where K is a constant. 

Using this result, the result that dE/d /jl is proportional to r 3 , and 
that v 2 is proportional to 1/r 3 , Equation (3A) is: 


y a r j 
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